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When both 5° and 5° can decay to the same final state composed 
>• ! of two vectors, the interference between them and those among three 

polarization states result in intricate phenomena. In this note we derive 



O ' the time and angular distributions for general B — > processes in 

2 ■ a form convenient for actual analyses. We then apply them to specific 

examples and clarify the CP violating parameters obtainable in the 
D*p and J/ipK* final states. The time distributions for the Z)*7r final 
><! ' states are also discussed. 

V 

§^ 1 Angular dependence 

The essential parts of this and next section can be found in many references [^. 
^ ; Here, we attempt to describe central concepts and derive critical expressions as 

;h ' simply as possible. 

1.1 Introduction 

We consider a two-body decay ^ a + 6 in the rest frame of the parent particle, 
where the spin state (J, M) of the parent particle and the heli cities of the 
daughters are given. The final state with a definite total angular momentum and 
definite helicities can be constructed as follows: In general, if \nX) is a state with 
total angular momentum A along the direction h, one can form a state with total 
angular momentum | J, M) where the quantization axis is taken as the z direction 
(i.e. in the lab frame), as 



|JM,A) = J dfiDfl^in) \nX) 
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with 

dn = d<P rfcos e , Dl^„, (n) = D^^^, (0, 0) (2) 

where {9, 0) is the polar coordinate of the direction n, and -D^,^,(n) is the rotation 
function, or the wave function of a top with total angular momentum | JM) and the 
component along n given by A which is also a good quantum number. 

Suppose IpAflAb) is the state in which particle a is moving in the p direction with 
helicity A^ and particle h is moving in the —p direction with helicity A;,: 

\p\ah) = \p\a)\~p\b). (3) 

Then, (|l]) with the identification 

n=p, \ = \a-\b, (4) 

gives the state with total angular momentum | JM) and total helicity A^ — A;, along 
the direction of a: 

\JM, = N j dpDi^,^^,^{p) \p\aK) , (5) 
where is a normalization factor. The ranges of the integration are 

-l<cos^<l, O<0<27r. (6) 
The possible values of the heclicities are constrained by 

|Aa-A,|<M, (7) 

which arises since the orbital angular momentum cannot have a component along 
the line of decay. The construction (j^) indicates that the amplitude for particle a 
to be in direction p is Df^ ^^^^^^^{p). 

Transformation of the state \JM, \a\b) under parity is given by 

P\JM, XaXb) = TraM-iy-'^-^'lJM, -Xa " A,) , (8) 

where Sa,b and T^a^b are the spins and intrinsic parities of the daughter particles, 
respectively. 

1.2 B ^ VaVi), helicity basis 

In B decays of the type B — > {V: a vector), such as 5+ '^K*+ and D*'^p+, 
we have 

J = M = 0, Sa = Sk = l, 7r„ = 7rfe = -l. (9) 
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The constraint with M = means Xa = Xb, and thus there are three possible 
hehcity states: 

(Aa,Afe) = (+l,+l), (0,0), or (-1,-1). (10) 
Accordingly, the final state can be written as 

\^j) = J2H,\h) (A = +1,0,-1), (11) 

A 

where is the amplitude for each helicity state, and we have defined 

^ |JM,+1 + 1), 
\fo) = \JM,00), {J = M = 0) (12) 

= |JM,-1-1). 

In terms of decay amplitude, one can write 

H^ = {fx\H,s\B), (13) 

where H^s is the effective Hamiltonian responsible for the decay. 
When the daughter particles subsquently decay as 

a ^ ai + a2 , 6 ^ 6i + 62 , (14) 

the construction (^ applies to each decay in its rest frame. The decay amplitude for 
ai to be in direction {9a, (f>a) in the rest frame of a and bi to be in direction {9^, (pb) 
in the rest frame of b is then (up to an overall constant) 

A = J2Hm ^rA.,-A., i<Pa, ^a, 0) (0^, ^b. 0) . (15) 

m 

The z axis in the rest frame of a is taken to be in the direction of p, and that in the 
rest frame of b is taken to be in the direction of —p; namely, each in the direction 
of the motion of the parent particle in the B frame. The definition of the azimuthal 
angles amounts to defining the phase convention for the helicity amplitudes H^- 
To be specific, we define that the x directions in the two frames are the same (see 
Figure Q). Using 

D^,^,(a,/5,7) = e-'^'^dl^^,iP)e-''-'^ , (16) 
the amplitude can be written as 

A = J:H^ e^'-^CA^.-A^, iOa)dl,^^_,^^ i9b) , (17) 

m 

with 

X = <i)a + (f)b (18) 

being the azimuthal angle from bi to ai measured counter-clock-wise looking down 
from the a side. 
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Figure 1: Definition of angles for B VoVh decay. 



1.3 Transversity basis 

Using the values @, the parity transformation (|]) reads 

P|/+i) = |/_i) P|/o) = |/o) P|/-l) = |/+l) 



(19) 



namely, the helicity-basis states | ± 1) are not parity eigenstates. However, we can 
construct parity eigenstates as 



and |/„) 



(20) 



The final state (|TTD can then be written as 

I%) = E^a|/a) 

A 



(A = I1,0,±) 



(21) 



with 



H,-H^ 



and Aq = Hq. 



(22) 



V2 



This basis is called the transversity basis P]. 

An often-used set of angles for the transversity basis can be obtained as follows: 
We note first that the angles (6'a, x) defined in the previous section is the polar 
coordinate of the Oi direction in the a rest frame where the 2;-direction is taken to 
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be opposite the direction of b in that frame and the x direction is taken to be in the 
decay plane of 6 — >• 6162 such that Px{bi) is positive. This defines a right-handed 
coordinate system where the y axis is perpendicular to the decay plane. We now 
define a new right handed system by 

x' = z, y' = X, z' = y, (23) 

where the z'-axis is now perpendicular to the b decay plane. Then, {9tr, Otr) is defined 
as the polar coordinate of the ai in this new system. Namely, {(ptr,4>tr) and {9a, x) 
are related by 

x' = sin 9tr cos (ptr = cos 9a = z 

y' = sin 9tr sin (ptr = sin 9a cos x = x . (24) 
z' = COS 9tr = sin 9a sin X = y 

These angles are shown for the case of D*p in Figure |^. Note, however, that one 
could also use the angles (x, 9a, 9b) for the transversity basis. 

1.3.1 B D*p+ (helicity) 

A full angular analysis of this mode has been presented at conferences [Q, but has 
not been published. Here, we consider the decay B D*p^ which is followed by 

D*^Dn, p+^7r+7r°. (25) 

We assign, 

a = D* , ai = D, aa = tt , b = p+ , 61 = 7r+ , 62 = 7r° . (26) 
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The decays of D* and p have only one hehcity state: 

^ai ~ = , Xb^ — Afej = . 

Thus, the general amplitude form (|r^) becomes 

m 

where we have relabeled the polar angles 

e = Oa (D*) , ^ = ef,{p). 

This can be rewritten as 

A = H+g+ + HoQo + H_ig_i 

where 

gj^i = -e*^ sin 9 sin ip 
go = cos 6* COS'?/' 

= -e~^^ sm9 sinip. 
2 

We have used 

together with 

The square of the amplitude is 

\Af = {J2Hmg^TiJ2Hngn) 

m n 

E\TT |2| |2 



Using the explicit forms for gm, the final distribution is 

r(x,^,^) = —\{\H+\'^ + \H_\^)sm'^9sin^ij + 4\Ho\\os^9cos^ip 
327r L 

+2{3?(i7+if*) cos2x - '^{H+Hl) sin2x} sin^^sin^ 

+ H^)H*) cos X - %{H+ - H^)H;) sin x} sin 29 sin 2^- 
where the normalization factor is chosen such that 

r dx t dcos9 dcostp T{x,0,'ip) = \H+\^ + + \Ho\^ . 

Jo J-i J~i 
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1.3.2 B D*p (transversity) 

Here, we can transform the D* side to transversity angles, or we could choose the 
p side. We arbitrarily choose D* side. We start from the amplitude (pOD and apply 



the transformations ( ^21) and (^J). We obtain 



^11 = -^(9+ + 9-) = ^ cos X sin 61 sin t/) = sin Otr sin sin i/j 

go = cos 6^ cos?/' = sin 9tr cos (ptr cos ip , (37) 

9± = ^(9+- 9-) = 73 sin X sin Osinip = cos Otr sin ^ 

to be used in 

A{(t)tr,Otr,i') =^A^g^{(t)tr,9tr,i') (m = ||,0, ±). (38) 

m 

Squaring this as before, the angular distribution becomes 

d^T{(ptr,Otr,i)) 9 / 2o • 2/) • 2 , • 2 / 

— — ; — ; = \A\\\ 2 sm. Sin (ptr Sin W 

dcptrd cos 6 trd cos ip 3277 V 

+ 1 y4_L 1^2 COS^ 6tr siv? l}) + \Aq\^A sir? 9tr COS^ (j)tr cos^ 

sin 6tr sin 2(j)tr sin 2iIj - V2%AqA±) sin 2etr cos (ptr sin 2ip 
-2'^{Ar,A^) sin 29^,. sin (ptr sin^ ip] (39) 



Integrating this over (ptr loses all interference effects among different polarization 
states: 



d''T{K.Otr.i^) ^ 9 
d COS 9 trd COS 16 



A\\\^ sin^ 9tr sin^ iIj+\Ao\^2 sin^ 9tr cos^ 'iIj+\A\\2 cos^ 9tr sin^ ip 

(40) 

At this point, we see that the even parity states (Ay and Aq) have sin^ 9tr distri- 
bution, and the odd parity state {A±) has cos^ 9tr distribution. Thus, plotting 9tr 
distribution only can separate even and odd parity components. On the other hand, 
both A\\ and A± are associated with sin^ipy ^ind thus ip distribution alone cannot 
separate different parity components. Further integrating over ip gives 

dT{(l)tr,9tr,-ip) 3 / 2 , I 4 |2\ • 2 /I , |/i2|o 2n \ /at\ 

= - (|A||| +\Ao\ )sin 9tr + \Aj^\2cos 9tr] ■ (41) 

d COS 9tr 4 V / 

1.3.3 B "^K* (helicity) 

The time-independent analysis has been performed by many experiments 0. We 
assign 

a = ^ , ai = i+ , a2 = i~ , b = K* , h = K , = vr . (42) 
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The decay K* Kit has only one hehcity state A^j 2 = 0- other hand, the 

final state of \E' — can have multiple helicity states because of the lepton spins. 
The actual helicity states, however, are restricted to only two due to the vector 
nature of the coupling that creates the lepton pair: 

(A,+ ,V) = ( + i,-i) or (43) 



2' 2' 
We have thus, 

Xa, - \a, = ±1 , K -\b,=0. (44) 

The final angular distribution is given by incoherent sum of the distributions for the 
two lepton helicity combinations: 

r(x,e,^) = |A(+if + |A(-i)r, (45) 

with 

A^'^=EHme''^^dUO)dU^), (46) 

m 

where we have relabeled the polar angles 

e = eam, i, = e,{K*). (47) 

The amplitude (^) has the form 

A^^^ =Y.Hmg^^\ (48) 

m 

where 

g^^i^ = — ^(1 + cos^^)e*^ sin-?/' ^f^^^ =— ^(1 — cos6')e*'^ sin 

5^0'*'^'* = sin 6^ cos , Qq =—-^sm6 cosip . (49) 

g^_^i^ = -^{1 — cos6)e~''^ sinip g^Si^ =^(1 + cos6')e~*'^ sin?/) 

The square of the amplitude (|48| ) is 

m n 

/ . l-'-'ml li/m I 

m 

+2 Y: {m*mHn) ^{9i^^*9i'^) - m*mHn) "^{gl^:^* gi'^)) • (50) 

Using the explicit forms for g!^\ the final distribution is 

T(x,9,i:) = — f(|iJ+|2 + |i7_n(l + cos^2)sin2V' + I^or4sin2^cos2^ 
647r L 

-2{^{HIH^) cos 2x + %HIH^) sin 2x} sin^ 6 sin^ tfj 

-{^{{H+ + H^yHo) COSX + '^{{H+ - H^yHo) sinx} sin 20 sin2V^] , (51) 



where the normahzation factor is chosen such that 



dx t dcosO r dcostp T{x,9,ij) = \H+\^ + \H_\'^ + \Ho\ 
Jo J-i J-i 



(52) 



1.3.4 B "^K* (transversity) 



One can use the relations (^2]) and (|2^) directly in the angular distribution ([5T|) to 
obtain 
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327r 



p(l — sin^ 9tr sin^ (f)tr) sin^ ip 



+ |/lo|^2(l — sin^ 9tr cos^ (j)tr) cos^ ip + sin 9^^, sinip'^ 

1 1 
— 3ft(y4j| y4o) sin 6^^^ sin 2(f)tr sin 2-?/^ + Q{AIA±)—^ sin 26*4^ cos sin 2ip 



V2 



V2 



+'^{A*^A±) sin 29tr sin 0^,, sin^ ip 
which is normalized as 

r2n 
/O 



(53) 



r27r rl rl 

d(f)tr J ^dcOs9tr J ^dcOS'^T{(f)tr,9tr,'^) = \Ai\\'^ +\A±\'^ +\Ao\'^ . (54) 



The transformation of angles can also be done at amplitude level. With the 
substitution of ampltudes (0), the amphtude for a given lepton total helicity A 
becomes 

(55) 



A^'^ =J2Amg^J:\ (m = |i,0,±) 



with 



9\ 



(+1) 



.+1) 
96 

(+1) 
9± 



— |(cosxcos^^ + isin x) sinip 
^ sin 9 cos Ip 

— I (cos X + ^ cos 9 sin x) sin ip 



9\ 



*• =|(cosxcos6' — isinx) sin?/) 

-1) 



^0 



(56) 



(57) 



--^ sin 9 cosip 

-|(cosx — 'icos6'sinx) sinip 

In order to apply the transformation from {9, x) to {9tr, (ptr), it is easier to multiply 
an overall phase factor which does not affect the final angular distribution. We take 

^(^_^^^ — cosx + ^ cos6' sin X 



sin 9tr 
cos X + ^ cos 9 sin x 



sin 9- 



for fli+i) 
9i 



9m 1 

for oi:^) . 



(58) 



tr 
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It is easy to see that these factors are indeed pure phases using the relations (^41): 



sin^ 9tr = sin^ 6tr cos^ (f)tr + sin^ 9tr sin^ 



Hr 



cos^ 9 sin^ 9 cos^ x 

1 •2, 2 a ■ 2 2, 2/)-2 

1 — sm X + COS 9 sm x = cos x + cos y sm x 



cos 6'^ + sin^ 9 cos^ x 
V ' 

(1 -cos2^)(l -sin^ x) 



— cos X + ^ cos sin x 



sin 6* 



COS X + i COS 6^ sin x 



sin 9, 



tr 



(59) 
(60) 



Muhiplying 



(+1) 



to 



(+1) 



we have 



9ll 9\\ ^ 



(+i)e(+i) a I ■ ■ \ ■ , -cosx + icos^sinx 

J — — f cos X cos + « sm X ) sm . (ol) 

2 sin 9tr 



Using 



(cos X cos 6* + i sin x) (— cos x + ^ cos 9 sin x) = ~ cos 6* — i sin 6*^ sin x cos x 



the phase-rotated is then 



^11 2 

Other functions are similarly obtained: 



— sin 9tr (cos (ptr + i COS 9tr sin ^t^) 



cos + i COS sin (j)tr) sin ^/^ . 



(62) 



(63) 



3\ 
9^0 



(+1) 



I (cos (ptr + i COS 9tT sin 0^^) sin ■?/' 

;^( — sin (ptr + ^ cos 9tr COS ^jr) COS 

(+1) I ■ a ■ I 
9± = ^ sm t^fr sm ip 



-1) 



=|(cos (j)tr — i COS ^^tr siu (j)ti) sin 

(70 =^{— SUKptr — i COS 9 tr COS (ptr) COS Ip 

g\_ =— I sin6'jr sin-?/^ 



(64) 



(65) 



These functions gives 

EAl^?i'¥ = |(l-sin2^i.sin2 



Hr, 



(A) 1 2 



'1 — sin^ 9tr cos^ (ptr) cos^ , X^A Wi. \ ' = \ sin^ 6'^,, sin^ ip 



siv? ip , 

(A) 1 2 



Ea = -i^ sin^ 9tr sin 20,, sin 2^^ , Ea = 



(A)* (A), 



(66) 



Ea Wr*9r) = , Ea ^(9^*9^) = ^ sin 29tr cos 0,, sin 2^ 
Ex^i9l'^*9^^) =0, EA^^(^?|f^*^7f ) = -isin2^,,sin0,,sin2 7/>, 
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which immediately leads to (p3D through (|50D where Hm are relaced by Am- Note 
that three of the combinations are zero; this arises from cancellations between the 
two lepton hehcities A = ±1. 

1.4 Charge conjugate decays 

For the charge conjugate decays {B decays), the rule for the definitions of angles 
is to start from the corresponding B decay, exchange paritcles and antiparticles, 
and then apply the definition of angles as if the daughter particles were the original 
particles from the B decay. For example, for the decays corresponding to assignment 
( P^ ) for B D*p'^, the particles in the decay B D*p~ are assined as 

a = D* , ai = D , 02 = vr , b = p~ , 61 = vr" , 62 = . (67) 

and the angles {6,x,'ip) are defined in the same way in terms of ai^2 and 61^2- In 
particular, the angle x is the azimuthal angle from bi to ai measured counter-clock- 
wise looking down from the a side. 

With this definition, the angular distribution is given by { \n\ ) with replacement 

Hx^Hx, (68) 

with 

Hx = {fx\H,s\B). (69) 
When CP is conserved in decay, then we can take (see Appendix) 

Hx = H_x (CP), (70) 

which holds to all orders in perturbation theory. In the literature, one sometimes 
encounters a CPT relation Hx = Hlx which is correct only to first order in pertur- 
bation theory. This CPT relation is thus not applicable to the decays of concern 
where the strong phases play inmportant role, since those phases are higher order 
effects. In terms of tranversity amplitudes, the CP relation (|68|) reads 

A||=A||, A^ = -A^, Ao = Ao. (CP) (71) 

Inspecting the expressions for the angular distribution, one notes that moving 
from B decay to B decay according to (^) or ( [7I| ) corresponds to changing x to — x 
for the helicity formulation, and 6tr —>■ n — Otr for the transversity formulation. These 
are nothing but the parity transformation (or equivalantly the mirror inversion) of 
the configuration. Namely, if one exchanges particles and antiparticles and take 
mirror inversion, then the resulting angular distribution is the correct one, which is 
to say that CP is conserved. 
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2 Time-dependence 



In this section, we will develop a formalism suited for neutral B decays to final 
states that are not CP eigenstates. In later sections, it will be applied to D*~^'k~ 
final state as well as each of the three polarization states of D*^p~ or "^K*. 

First, let us recall the time evolution of pure and B^ states. Assuming CPT, 
the physical states Ba and B}, can be written as 

Ba=pB^ + qB° (m„,7,) 

B,^pB^-qB^ K,7(.) ' ^ ' 

where ^ and 7^,5 are the masses and decay rates of the corresponding physical 
states. Theoretically and experimentally, \p\ = \q\ within error of order 1%. Here, we 
assume \p\ = \q\ which makes p/q a. pure phase factor. The lowest order estimation 
gives (see Appendix) 

(73) 



which corresponds to the choice of the CP phase of the neutral B meson given by 
CP\B^) =r]B\B^) , CP\B^) =7]%\B'^) , with rjB^l. (74) 
The above value of p/q is for the case Ba is heavier than Bi,: 

ma > nib. (75) 

The physical states evolve as 

Ba S„e-^"^«*-^* , Bb Bbe-'""''-^' . (76) 

Hereafter, we will assume that the decay rates of the two physical states are the 
same 

7a = 7b = 7 • (77) 

then, the factor e~^* decouples from all amplitudes, which we will drop for now and 
restore it at the end. We also separate an overall phase factor exp(— i "^''^"^'' t) and 
discard it since such overall phase factors do not affect measurable quantities. Then 
the evolutions of Ba^b can be simplified as 

Ba^Bae-^'f^ Bb^BbC^'f' (xe-i*) , (78) 

with 

5m = ma — mb > . (79) 
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Then, the time evolutions of pure and can be obtained by solving (|72D for 
B^ and B'^ and then applying the time evolutions above: 

B^ = — ^ —{ BaC'—^ + Bhe'~^) = B° COS -5°zsin . (80) 

pB^ + qB^ pB^-qB^ 

The time evolution of B is similarly obtained. Restoring the decay factor e~2 , 

dO -it /^dO '^"^^ ^ dO- • "^"^A 
^ e 2^ \ B cos -B z sm , 

5° e-i* (5° cos ^ - ^5°. sin ^) . (81) 

We now consider the decay amplitudes for a pure B^ or state at t = to 
decay to a final state / or its charge conjugate state / at time t. The final state 
could be D^*^^n^ or any given polarization state of D*^ or "^K*'^. Define four 
instantaneous dedcay amplitudes by 

a = Amp{B° f) 
a = Amp^B"^ /) 

h = Amp{B'' f) ' ^ > 

b = AmpiB^ f) 

For / = D^*^~n~^, for example, a and a are the favored amplitudes and b and b are 
the suppressed amplitudes. Then, (|8l|) gives 



^3) 



Aso^fit) = e-¥ (^acos^ - |Hsin^)= e-^^a (cos^ -pisin^ 
AB^^fit) = e-^* (acos^-26isin^)=e-^*a(cos^-p2sin^ 
ABO^jit) = e-2* (^6cos^ - %isin^)= e-i*a(p cos ^ -2sin^ 
Aso^jit) = e-2* (^6cos^ - Raisin e-i*a (pcos^-zsin^ 

with 

p = — , p = — . 84 

pa qa 

For the bottom two amplitudes (the 'suppressed' decays), we have ignored overall 
phase factors p/q and q/p for the second equalities. 

At this point, we can see the relation between the 'suppressed' and 'favored' 
modes; namely, up to an overall phase, 6mt Smt + vr transforms A^o^flt) to 
ABO^fit) and ABO^jit) to ABO^f(t). Equivalently, in the expressions of decay rates, 

(cos 6mt, sin drnt) ^ (— cos 6mt, — sin 6mt) (85) 
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transforms between a suppressed mode and its favoed mode with the same final 
state. Also, p/q is the complex conjugate of q/p (within the approximation that 
\p\ = \q\), and as we will see more explicity later, the weak phase of b/a is the 
complex conjugate of that of b/a with the rest being the 'strong phase' which is 
common to both. Thus, 

(weak phase) —(weak phase) (86) 

keeping the strong phase the same transforms between a decay and the corre- 
sponding decay (both 'suppressed' or both 'favored') apart from the difference 
between a and a. Often \a\ and \d\ are the same and if so the above transformation 
is exact in the decay rates. When we extend the above time-dependent amplitudes 
to include interferences between polarizations, the rule between the same final state 
(|85|) still holds, but the relation between 5° and 5° ( ^61) does not hold in the helicity 
basis. We will see, however, that it holds in the transversity basis. 
The time dependent rates are obtained by squaring (pS]): 

[(1 + IpP) + (1 - IpP) cos6mt + 2Qp sin^mt] 
a|2^ [(1 + |p|2) + (1 - |p|2) cos6mt + 25>p sindmt] 

(87) 

a|2^ [(1 + IpH - (1 - |p|2) cos6mt - 2Qp sm6mt] 
ap^^ [(1 + IpP) — (1 — IpP) cos6mt — 2Q=p sin^mt] 



In deriving this formula, we have assumed CPT in the mixing and that 7a = 'fb- 
Otherwise, it is general; in particular, there could be direct CP violations in any of 
the decay amplitudes such as \a\ ^ \a\ etc. 

On T(4S'), one would fiavor-tag the other side by, say, a lepton. If the tag side 
decays to i~ at proper time Uag, the quantum correlation is such that the signal 
side is pure 5° at the same proper time tgig = Uag and proceed to evolve as usual 
from that time on. Thus, for tgig > Uag, the decay distribution is given simply by 
the replacement 

t — > At = tsig — ttag ■ (88) 

For tsig < Uag, all that is needed is to put absolute value on At of the decay factor 
g-7Ai ]Sfa]2iely, (|87|) becomes the distributions on T(45') with the replacement 

jt 7| At| and Smt 6mAt . (89) 
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Figure 3: Diagrams for D n 



Explicitly, 
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pp) cos^mAt + sinSmAt] 
pp) cosSmAt + 253p sin5mAt] 

, (90) 

p| ) cos5mAt — 2Qp sm6mAt] 
pp) cosSmAt — 253p sin5mAt] 



where F^- j(At) denotes the decay rate for one side decaying to a final state / while 
the opposite side is tagged by a negative lepton (or tagged as by any other 
method), etc. 

2.1 L»W-7r+ 

Earlier studies of this mode can be found in Ref. 0. Diagrams for B^ —>■ D^ir^ are 
shown in Figure ^. In addition to dominant tree diagrams, annihilation diagrams 
may have non-negligible contribution. Also, there may be final-state rescattering 
Di*)On^ L)(*)-7r+. The CKM factor of these processes, however, is the same as 
that of the corresponding tree diagram for the same final state, and thus it does not 
affect the following formulation. Penguins should result in even number of charms; 
thus, penguins do not contribute. 

With the definitions / = D^n^ and / = D^n^ , the four amplitudes of (|8^) can 
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be written as 



a = 

a = 

h = 

h = 



Amp{B° - 


D- 


-7r+ 


) = KFc 


Amp{B'^ - 






) = KFc 


Amp{B'^ - 






) = KFu 


Amp{B^ - 


D- 


-7r+ 


) = KFu 



with 



K 



Vctv:, 



(91) 



where we have separated the CKM factors K*}^ and called the rest Fc^u which in- 
clude strong phases as well as decay constants and form factors (if factorization is 
assumed). We assume that the CP violation is solely through the weak phases that 
appear in (^); as a consequence we can show that (see Appendix) 



We then have 
Using (^) 



\a\ = \a\ , \b\ = \b\ . 
and (|8^ as well as (|9^), the value of p defined in (|8^ is then 

pb pKFu 



(92) 
(93) 



P 



VtWtbVub^cd Fu ^ 



qa q X^F^ 

where we have defined r = |p| and (pp = argp. With the definitions of and 03 



(94) 



)i = arg 



h = arg 



-VcdV:, 



we have 
arg 



\ vuv:,v,,v:j 

Then, p can be written as 
with 

Similarly, one obtains 
Note that we have \p\ = \p 

Vub^cdFu 

VcbV:, 



arg 



VcdV:, v,,vx VudV:, 



-VtdVg -VtdVl -VrdVX 



+ 



re 



^3 , 



tdVtb-VcdVcb, 



= arg — . 



p = re-'^'^^-^^ . 
= r. The value of r is roughly 

OAX^ ~ 0.02 . (A ~ 0.22 : Cabibbo factor) 



(95) 

(96) 

(97) 
(98) 
(99) 

(100) 
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With '^p = — rsin(0tt, — 6) and = rsin(0„, + 6), the four decay rates (p^) 
becomes 

^BO^D-TT+ii) =|ap^^ [(1 + '"^) + (1 ~ '"^) COS 6mt — 2rsm{(f)w — 6) sm6mt] 

^BO^D+TT-it) =|ap^ [(1 + r^) + (1 - r^) cos 6mt + 2rsm(0^ + 5) sin^mt] 

rB0^D+7r-(^) =kP^-Y~ [(1 + ~ (1 ~ cos6mt — 2rsin(0^ + 6) sm6mt] 

T§o^n-K+(t) =kP^-Y- [(1 + r"^) - (1 - r'^) cos6mt + 2rsm{(f)y, - 6) sm6mt] 

(101) 

where have used \a\ = \a\. Note that T§o^£)~T,+ {t) (suppressed) is obtained from 
TBO^D-Tr+(t) (favored) and Fbo^d+vt- (t) (suppressed) is obtained from r5o^£)+^- (t) 
(favored) by the transformation (pSj), and within the two suppressed modes and 
within the favored modes, the expresssions are related by ( ^61) namely 0^ ^ —(pw 
The CP violating parameters that can be extracted from these distributions are 

rsin(0^ — 5) and rsin(0^ + 5). (102) 

Note that the two extractable paramters are always multiplied with r, and the 
value of r cannot be obtained by the fit. As discussed earlier, the corresponding 
distributions on T{4S) are obtained by replacements 7^ — > 7|At| and 6mt — > 5m At. 
The first paramter rsin(</)^ — 6) can be obtained through the asymmetry between 
positive and negative At of F^- £)-^+(At) (favored) or Ti+^D-T,+ {At) (suppressed), 
and the second paramter rsin(0^ + 6) is similarly obtained through £)+^-(t) 
(favored) or T^- £)+^-(t) (suppressed). This feature that single mode can give a CP 
violating parameter through asymmetry between positive and negative At is unique 
to T(4S'). In fact, most of the information on CP violation is in such asymmetries. 
If we define 

STxi\At\) = Tx{At) - Txi-At) , (103) 

we have 

-6Te-^D-^+{\At\) =6Te+^D-^+{\At\) = Nr sm{(f)^ - S)e-^' sm{Smt) , (104) 
^r^+,D+^-(|At|) = -6Te-,D+n-i\At\) = Nr sm{(f)^ + 5)e-^' sm{6mt) , (105) 

where is a common normalization factor which is known. 

Now we derive the corresponding time-integrated expressions. We use following 
integrals. 

lo^e-^'dt =i, (106) 

/n°° e-T* sin 5mt dt = -— ^ , (107) 

r e-^* COS 6mt dt = . (108) 

7 1 + x"' 
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Here x = Sm/'y. The time-integrated decay rates become 



r(5° - 


D 




l^|2 


[(l + r2) + 


1 
i 






2rx 


sin(( 




27 


1 


+ x"^ 




1 + x^ 








|a|2 


[(l + r2) + 


1 




+ 


2rx 


sin(( 


+ 5)] 




1 


+ 


1 + x2 


r(5° - 


^D' 




|a|2 




1 






2rx 


sin(( 


K + 5)] 


" 27 


1 


+ 




1 + x2 


r(5° - 


D- 




|a|2 


+ 


1 




+ 


2rx 


sin(( 






1 


+ X2 


1 +X2 



If we set (5 = for simplicity, we see that the information on sin(0^) is in the 
asymmetry between the top two rates (the favored modes) or in the asymmetry 
between the bottom two rates (the suppressed modes). The absolute amount of the 
differnce is the same for both cases, but the total rate is about 5 times larger for the 
favored modes compared to the suppressed modes. It means that the significance 
(number of sigmas) is times smaller for the favored modes. Thus, most of the 
information is contained in the suppressed modes. 

The expressions ( |101| ) and ( |109|) are valid also for / = D*^7r+, D^p^. When 
there are more than one polarization states as in D*~ p^, there is extra effect due to 
interferences between different polarization states, which we will discuss next. 



2.2 D'-p+ 

This mode was first stdied in detail in Ref. |0. We first note that the expressions for 
the time dependent amplitudes ( [53| ) are still valid when applied to each polarization 
state: 

ABO^f^ (t) =e-i* (oa cos ^ - ^hi sin ^) = e'i* ax (cos - pxi sin ^) 

Aso^f^it) =e-i* (ciA cos ^ - J6a« sin ^) = e'^* cia (cos ^ - pa « sin ^ 

Abo^j^ {t) =e-i* (6a cos ^ - ^dxi sin ^) = e'^* ax {px cos ^ - i sin ^ 

Aso^j^it) =e-i* (6a cos ^ - 2aAi sin ^) = e'i* oa (pa cos ^ - i sin ^) 

(110) 
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where 





fx] 


) = KFcX 


dx= Amp{B^ - 


fx, 


) = KFcx 


h\= Amp{B^ - 


fx, 


) = Ki^uX 


bx= AmplB^ - 


fx] 


) = ^uFuX 


qbx 

P\ - 

pax 


Px 


_ pbx 
~ qax' 



(111) 



and 

(112) 

Each of ( |110| ) gives the polarization amphtudes to a given final state at time t. 
Then, the angular distribution of pure i?'' at t = decaying to / at time t is simply 
obtained by replacing Hx or Ax by ABO^f^{t) in ( P5| ) or (P5|). For B^{t = 0) ^ 



/, for example, the time-dependent angular distribution is given by (|39|) with the 
replacement 

\Ax\' ^ \As.^fM\ 
^{AlA,)^^{A*j,._^f^^{t)As^^fM). 

»(A*/l^) -.$5(A^„^^„(t)A^o^^Jt)), 

^{AlA^) ^{Al,^f^^{t)Aso^f^{t)) . (113) 
Or the decay amplitudes are obtained from (^) by the same replacement: 

ABO^f{il,t) = 2^e 2 ax y^os— pA^sm^— j gx{^l) 

ABO^f{U,t) = 2^e 2 ax y^os— pA^sm^— j gx{^l) 

A .\ ^ ^:Lt- (- . . 6mt\ 

ABO^f{ll,t) = l^e 2 axypx cos ^ ism^- j gx{il) 

A^o^fin, t) =J2 e"^* ax (^px cos - i sin gxi^) , (114) 

where Q = {(f)tr, Otr, i^) or (x, 6, ip), and gx's are given by (^) or (|3T|) . Here, the final 
states are / = D*^ and / = D*^p^ . Note that one could use angles {(j)tr, Otr, ip) or 
(x, 0, ip) for the tranversity amplitudes (for that matter, for the helicity amplitudes 
also - we just have not provided g{4>tr, (^tr, i') for the helicity amplitudes). 

Since Fux and Fcx are nothing but the polarization amplitudes apart from the 
CP violating phases, they themselves should satisfy the CP relations (^) and (fn]) 
(see Appendix): 

Fgx = F,^x (helicity), (115) 
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Fq\\ = , Fqo = FgQ , Fg_i_ = -Fu± , (tianversity) (116) 
where q = u or c. 

2.2.1 Helicity basis 



With the relations (|115|) for hehcity basis, the decay amphtudes can be written as 



(hehcity) (117) 

bx= A*FuA , = ^uFuX ■ 

This gives 

\ax\ = \a^x\ , \bx\ = \b-x\ ■ (118) 
And the parameters px and px becomes 

Pa = - , p-x = - . • (119) 

P K^cX Q -^c-TcX 

Then, the same procedure that led to (p7| ) and ( p9D allows one to write 

= rxe'^'^^+^^^ , p_A = rxe-'^'^^'^^^ , (120) 

where 

IpaI = |P-a| = rx , (121) 

and 

Sx^arg^. (122) 

2.2.2 Trans versity basis 



Using the relations ( 116 ) for transversity, we can write 



ax— A*FcA , a A — ^aAc-^cA , 

(transversity) (123) 

bx= A*F„A , bx = ^.x^uFux , 

where 

. _f 1 (A = 11,0) 

^^-\-l (A = ±) • ^^^^^ 
We will use the tranveristy basis for the rest of this section. Clearly, we have 

\ax\ = \ax\ and |6a| = I^aI , (125) 
20 



and the procedure semilar to that led to ( pTl) and (|99|) gives 



Pa = i\r\e 



Pa = ^'^Ae 



where 
and 



IPaI 



IPaI 



6\ = are 



uX 



(126) 

(127) 
(128) 



cA 



Let's evaluate the expilict decay rates; namely, the coefficients given in ( p.l3| ). 
Note that px and px are related by (p^j ^ —<t>w Togethter with (|85|) , all we need is to 
evaluate one of the four modes which we take to be the favored mode — »• D*~p^. 
Calculation is straightforward and we obtain (apart from the common factor e~^*/2) 



+ 



\ax\ [l + r^) + [1 — r^j cosomt 

-2^xrx sin(0^ - 6x) sin Smt , 
3?(ai|ao)(l + r||rocos(5|| - 5o)) + 0^(011 ao)?'Fo sin(5|| -60) 
3?(ajJao)(l — r||ro cos(5|| — Sq)) — 53(a|Jao)rj|ro sin((5|| —60] 
5R(ai|ao)(r|| sin(0^ - (5||) + tq sin(0^ - 5o)) 
+$5(ai|ao)(rj| cos(0^ — (5||) — vq cos{(j)w — 60)) sin Smt , 
53(a*ax)(l — rer± cos((5e — S±)) + ^{ala±)rer± sm{6e — Sj 
53(a*a_L)(l + rer_L cos(5e — (^x)) — ^{ala±)rer± sm{6e — 5j 
5>(a*a_L)(^eSin(0^ - (5e) - r_L sin(0^ - 6±)) 



cos 5mt 



(129) 



cos 6mt 



(a*a_L)(re cos(0^ - 6e) + rj_ cos(0^ - 6j_)) 



sin 5mt , 



where A = (||,0, ±), e = (||,0), and the suppressed modes for the same final states 
are obtained by the transformation 6mt 6mt + tt or (P5|), and among the two 
suppressed or among the two favored modes, the decay and the decay are 



related by 



—611,. The distribution 



with these replacements then gives the 



desired time-dependent angular distributions. 



2.2.3 Fit parameters 

Squares of the amplitudes (|114|) give the rates, and with complex functions in pro- 
graming language, these expressions are all needed to perform the fit. The fit pa- 
rameters are ax, dx, Px, and px- Note that only the relative phases matter among 
ax and among dx', namely, one can set oq = real and do = real, for example. In 
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addition, IaaI = |o-a| (helicity) or \a\\ = \ax\ (transversity) reduces the number of 
degrees of freedom by 3 in each basis. Furthermore, there are phase relations in 

^ = ^ (A = ±l), or ^ = 6^ (A=||,±), (130) 

which reduces 2 degrees of freedom. Thus, there are 5 degrees of freedom in ax and 
dx including the overall normalizations. One may parametrize, for example, as 



a+i 




a+i e*"^ 


'+1 


a+i 




a_i e 






l^ol 




ao 




|ao| 


a_i 




a_i e*"^ 


'-1 


a_i 




|a+i|e 


a\\ = 


\a 








a\\ \e 




ao = 


«o 


, ao 




ao 





(helicity) , (131) 



(transversity) . (132) 



a± = |a_L|e*'^-L a_L = -|a_L|e''P^ 



Also px and px are constrained by the expression ( |112|) ; namely, we actually fit rx, 
Sx, and 0^, which amounts to 7 degrees of freedom. The total number of degrees 
freedom is thus 5 + 7 = 12 including the overall normalization. 



2.3 Time dependent angular distribution for 

The only one relevant final state to be considered is where K*'^ decays to 

Ksn'^. We denote the final state as fx = {'ifKg^)x, where A could be for hehcity 
basis or transveristy basis. The particle assignments are 

a = ^ , ai = , as = r , b = K*\ h = Ks , aa = 7r° . (133) 

All we need is the amplitudes for each polarization (helicty basis or transversity 
basis) at time t when the B meson was pure 5° or 5° at t = 0. Then, we can 
use the distributions (0) and ([53|) to obtain the angular distribution at that time. 
Incoherent sum over the two possible helicity states of the \E' decay is already taken 
into account in those angular distributions. 

Since we are dealing with only one final state (apart from polarization), the first 
and the last of (^) will do: 

Aso^f^ (t) =e-i*aA (cos - px^ sin ^) 

' ' , (134) 

A Bo^f, (t) =e- ^ *aA (pa cos ^ - z sin ^) 

where 

ax = AmpiB'^{^Kf)x), hx = Amp{B^ Kf)x) . Px = ^ ■ (135) 

pax 
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The Pa parameter is then 

(136) 



q{Ks 


((^i?*0)A 


Hcs 


BO) 









Eq. ( [73D gives g/p, and using (|185|) of Appendix, 

(137) 



can be obtained as in the case of B where we have ignored the small deviation of 
\<1k/pk\ from unity. Assuming that the color-suppressed tree diagram dominates 
the amplitudes a and 6, or assumming that penguin and other contributions do not 
modify the weak phase significantly, 

((^ir*0)A|ifcff|50) = KtKsi^A , ((^ir*°)A|i^cfr|5°) = KfeKt^A . (138) 

By the similar argument that led to the CP relations ( |115D and ( |116| ), Fa and Fa 
are related by 

FA = i^-A (helicity), (139) 

F|l = i^i , Fo = i^o , F^ = -F± , (tranversity) (140) 

Let's use the transversity basis for the rest of this section. Then, the amplitudes 
given by ( |138| ) togehter with the CP relation above gives 

{(9K-WH,„\B") V^V:, 
where ^a is the sign defined by ( |124|) . Combining all ingredients, p\ becomes 

- rvw;J [v^] K^'w:.) - l^v^j / ■ 



With the definition of (pi (0), we can write 

= -^^e'^'^^ (transversity) . (143) 

Recall that the value of p for the gold-plated final state was e"^**^^; namely, 

the transverse polarization has the same time-dependent CP asymmetry as the 
final state, and A\\ and Aq states have the CP asymmetry opposite to that 
of '^Ks- These arguments are valid when a given polarization state dominates the 
final state and when integrated over the angular distribution. 
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The angular distribution is given by the expression ( ^3]) with the coefficient 
replaced according to (|113|) . Explicitly, 



I^Ap |aAp(l ± sin20i sin 5mt) 
3fJ(v4| v4o) ^ 3?(a|Jao) (1 ± sin 20i sin Smt) (144) 
^{AIA±) ±53(a*a_L) cos(5mt =1= 3?(a*a|j) cos20i sin(5mt 



where the upper sign is for 5*^ "^Kg^, the bottom sign is for "^K^ 



S 5 

A = (II, 0, _L) and e stands for || or 0. They are related by the transformation ( p5D as 
expected. Note that cos20i can be obtained by these angular distributions, which 
helps to resolve the discrete ambiguity of 0i. 

3 Appendix 

3.1 CP relations 



We will hereby derive the CP relations (|70|) and (|115|) . Suppose the effective Hamil- 
tonian commutes with CP: 

IH,sP = H,s, (145) 

where 

I =CP. (146) 

For example, the Hamiltonian for the tree diagram of D^n^ is (up to a 

constant) Q 

^eff = KK + (h.c.) (147) 

with 

hc = j'^dt j d^x{ch)^{duY , (148) 



rT 

-T 

where {qq')^ is a color-singlet V — A current which is a function of space-time: 



{qq'ix)), ^ rixh.il - l5)q''ix) . (x = (t, x)) (149) 



with a being the color index. The CKM factor Ac is defined in (|8^). The CP phases 
of quark fields are taken as 

r/, = l, (150) 

where the CP phase is defined by 

^kpa) = ^?g|g-p,a) , (151) 



^What we are calling H^s is actually the S operator. 
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where p is the momentum and a is the spin component along z. With this choice of 
CP phase, one can show that (see, for example, Ref |^) 

J(gg'(x))^/t = -{qq\x')Y^ = -(g'g(x'))'^ (x' = (t, -x)) . (152) 

Note that the Lorentz index /x changed from subscript to superscript. We then have 
after the integration over space 



£x{chY{du)f, 



(153) 



which leads to 

IhJ^ = h\. (154) 

Similarly, we can show 

IhlP = h,. (155) 

This makes (|145|) hold if Ac is real. In general, ifgfr includes strong interaction that 
results in phase shifts. Still, it can be writen in the form ( |147| ) and that it would 
be invariant under CP if the CKM factors are real; namely, ( |154| ) and ( [155| ) are 
satisfied. 

The helicity states of i? — >■ a + h transforms under CP as (see, for example, 
Ref §) 

/ I JM, ah) = r]aVb{-y-'''-'VM, -A, - A^; db) , (156) 

where J = and = 1 for our case, and rja and rjh are the CP phases of a and b 
respectively: 

C P\a; p, a) = r]a\a; -p, a) , C P\b; p, a) = r]b\b; -p, a) , (157) 

where a is the z-component of spin. If a or 6 are not self-conjugate, then their CP 
phases are cancelled when the value of p is calculated or relation between px and px 
is evaluated. For a self-conjuagte particle, the CP phase does matter. However, CP 
of relvant spin-1 particles, such as any known spin-1 (cc) states, p°, a°, u, 0, etc. 
are all +1. Thus, we take 77^776 to be +1 keeping in mind that if any of the spin-1 
particles are self-conjugate and CP— then it has to be included in the sign. Thus, 
in terms of our short notation, the above relation becomes 

/|/a) = I/-a). (158) 

Then, the helicity amplitude transforms as (with 77^ = 1) 

Hx = {fx\H,s\B') (159) 
= ifxlI^IH^ ^B^=H^x, (160) 
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which proves ([70D. 

The proof of the relation (|115|) starts from reahzing that the relevant effective 
Hamiltonian can be written as 

H,s = (Xchc + \uK) + {Khl + Khi) , (161) 

where the second part is just the h.c. of the first part to make the whole Hermitian, 
and Ac = VcbV*^ and Xu = VubV*^ as before. The term Xchc includes ab ^ c transition 
and creation of a ud pair, and X*hl includes ab-^c transition and creation of a du 
pair, etc. Again, contain the effect of strong interaction to all order. 

Assuming that CP violation occurs solely through the complex CKM phases, 
we should have 



IhlP = h,, IhiP = K, (162) 

which makes Hes invariant under CP if Ac,n are real. Then, ax defined in ( |111| ) with 
/ = D*~ can be written as 

ax ^ {fx\H,K\B') = {fx\X:hl\B') = X:{fx\hl\B') . (163) 

Comparing with (|1 1 1|) indentifies Fcx as 

F^x = ifMlB') . (164) 

Similarly, 

F,x={fx\hc\B''). (165) 
Then, we have (again with 77^ = 1) 

F,x = {f>^lllJ^r\B^=F,^x. {QED) (166) 

The proof of F^x = F^-x (helicity) proceeds the same way. 

The relations (|92D is proved similarly. Here, the effective hamiltonian is again 
written in the form ( |161D and is invariant under CP if the CKM factors are real. 
With J = si = S2 = A = in ( |156|) , the transformation of the final state / = 
7r+ under CP is 

I\f) = VDriAf). I\f)=V*DV:\f), (167) 

with proper choice of CP phases (you can set A = in ( |158| )). The quantities Fc,„ 
and Fc^u sue identified as 

F, = {f\hl\B^), F, = {f\K\B'), 

F^ = {f\K\B'), F^ = {f\hi\B'). (168) 
The final result is obtained by simply setting A = in (|166|) : 

Fc = Fc, and similarly Fu = F^. (169) 
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3.2 Derivation of 

The two mass-eigenstates are the eigenvectors of the Schrodinger equation in the 
B^-B^ space: 

= , (170) 

at 

where ^ is a two component vector and if is a 2 x 2 matrix in the B^-B'^ space: 



where M and T are hcrmitian matrixes. When 7^ = 7^, the decay part decouples 
and the mass matrix part (mixing part) only should be considered; thus, we will 
drop r. The CPT invariance allows one to write 

«=^=(r. «)■ (^-'^)- (1^2) 

In particular, 

II = {B'\H,s\B^) , li* = {B'\H,s\B') . (173) 

The eigenvalues are 

^''("m^' ™-a)="' ^ A = m±|H. (174) 

Let's define the eigenvector for the heavier of the two to be pB^ + qB^, which then 
should satisfy 

:)(:)=('"+i''i)(:)- (^^^) 

The top component (the B^ coefficient) of this equation gives 

mp + iiq — mp + — > - = . (176) 

Q W 

On the other hand, the B° B° transition is caused by the box diagram at the 
lowest order whose effective Hamiltonian can be written as 

H,fi = iVt,V:,nbd^S + {h-c.), (177) 

where hf^^^^i is the effective hamitonian that transforms B^ to B^ and itself trans- 
forms under CP as 

I hbd-^db = ^Id^db ' ^ ^Id^db ^ ; (178) 



r 
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namely, Hes is invariant under CP if it were not for the CKM phases. Then the 
off-diagonal elements of H are related by CP as (with 77^ = 1) 



*\2 

td) 



{B'\H,s\B') 



*\2 

td) 



^* . (179) 



Thus, 



td) 



{v:,v,,Y i/i| {v:,vu) 

Using ( |1 76|) , we see that pB^ + qB^ with 



;i80) 



:i8ii 



td) 



This method is simple and elegant but cannot define the sign; in order to do so, one 
needs to actually evaluate /i [ p^ : 

^ = -^i^f^msBBV2So{V,lV,,r (182) 

where fs is the decay constant of 772 > is a QCD correction factor, 6*0 > is 
a function of the top quark mass, and Bb is the 'bag factor' of the B meson which 
is believed to be positive. Then, p/q is now 



P_ VtbVtd 

q VtlVtd 



:i83) 



and these p and q makes + the heavier of the two mass eigenstates. 

In the neutral K system, we define px and qx in parallel to the B system; 
namely, the heavier (Kl) is defined to be PrK^ + qxK^- Thus, Ks is 

Ks = PKK''-qKK'. (184) 
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Even though there is some comphcation due to the hfetime difference; the situation 
for the phase of Pk/qk is essentially the same and to a good accuracy it is given by 
applying t ^ c and 6 — > s to (|183| ): 



cd 



:i85) 



The conventions used here for p and q are not the same as those in Ref. |11 . 
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